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ABSTRACT 

The rotation rates of six weakly-magnetic neutron stars accreting in low-mass X-ray binaries have most likely 
been measured by Type I X-ray burst observations with the Rossi X-Ray Timing Explorer Proportional Counter 
Array. The phenomenology of the nearly coherent oscillations detected during the few seconds of thermonu- 
clear burning is most simply understood as rotational modulation of brightness asymmetries on the neutron star 
surface. We show that, as suggested by Strohmayer and colleagues, the frequency changes of 1-2 Hz observed 
during bursts are consistent with angular momentum conservation as the burning shell hydrostatically expands 
and contracts during the burst. We calculate how vertical heat propagation through the radiative outer layers of 
the atmosphere and convection affect the coherence of the oscillation. We show that the evolution and coher- 
ence of the rotational profile depends strongly on whether the burning layers are composed of pure helium or 
mixed hydrogen/helium. Our results help explain the absence (presence) of oscillations from hydrogen-burning 
(helium-rich) bursts that was found by Muno and collaborators. 

We also investigate angular momentum transport within the burning layers and address the recoupling of the 
burning layers with the star We show that the Kelvin-Helmholtz instability is quenched by the strong stratification, 
and that mixing between the burning fuel and underlying ashes by the baroclinic instability does not occur. How- 
ever, the baroclinic instability may have time to operate within the differentially rotating burning layer, potentially 
bringing it into rigid rotation. 



Subject headings: accretion, accretion disks ■ 
bursts 

1. INTRODUCTION 



nuclear reactions — stars: neutron — stars: rotation — X-rays: 



Type I X-ray bursts have been long understood as thermonu- 
clear flashes on the surfaces of neutron stars accreting at rates 
of 10""Mo yr"' <M < \Q~^Mq yr"' in low mass X-ray bina- 
ries (LMXBs). The accreted hydrogen and helium accumulates 
on the surface of the neutron star and periodically ignites and 
burns. Thermonuclear flash models successfully explain the 
burst recurrence times (hours to days), energetics (^ lO^'^ergs), 
and durations (^ 10-lOOs) (Lewin, van Paradijs, & Taam 1993; 
Bildsten 1998), though many quantitative comparisons to ob- 
servations are less successful (for example, see Fujimoto et al. 
1987; Bildsten 2000). 

Evolutionary scenarios connecting the neutron stars in 
LMXBs to the millisecond radio pulsars (see Bhattacharya 
1995 for a review) predict that neutron stars in LMXBs should 
be spinning rapidly. This has been confirmed for one sys- 
tem, the V, = 401 Hz accreting pulsar SAX J1808.4-3658 (Wi- 
jnands & van der KHs 1998; Chaki-abaity & Morgan 1998), in 
which the neutron star magnetic field {B ^ lO'^-lO'' G; Psaltis 
& Chakrabarty 1999) channels the accretion flow onto the mag- 
netic polar caps, creating an asymmetry which is modulated by 
rotation. However, most neutron stars in LMXB's show no evi- 
dence for coherent periodicity in the persistent emission, imply- 
ing that the neutron stars do not possess magnetic fields strong 
enough to make a permanent asymmetry (B <C 10'"G). 

Type I X-ray bursts have provided a new way to determine 
the spin of these neutron stars. Observations with the Propor- 



tional Counter Array (PCA) on the Rossi X-Ray Timing Ex- 
plorer (RXTE) of neutron stars in six LMXBs have shown co- 
herent oscillations during Type I X-ray bursts, with frequencies 
1^0 that range from 300 to 600 Hz (see Table |l]). The simplest 
interpretation is that the burning is not spherically symmetric, 
providing a temporary asymmetry on the neutron star that al- 
lows for a direct measurement of rotation. The coherent nature 
of the periodicities (Q ^ 300), large modulation amplitudes and 
stability of the frequency over at least a year support its inter- 
pretation as the neutron star spin (Strohmayer 1999b, and ref- 
erences therein). 

It is expected theoretically that the burning should not be 
spherically symmetric. Joss (1978) and Shara (1982) suggested 
that ignition of a burst occurs at a local spot on the star, and not 
simultaneously over the whole surface. This is because it takes 
hours to days to accumulate the fuel, but only a few seconds 
for the thermal instability to grow. Simultaneous ignition thus 
requires synchronization of the thermal state of the accreted en- 
velope to one part in lO^-lO'* over the surface. More likely 
is that ignition is local and a burning front then spreads later- 
ally (Fryxell & Woosley 1982; Nozakura, Ikeuchi, & Fujimoto 
1984; Bildsten 1995), burning the rest of the accreted fuel and 
creating a temporary brightness asymmetry on the neutron star 
surface (Schoelkopf & Kefley 1991; Bildsten 1995). Consistent 
with the picture of a spreading burning front, the pulsation am- 
plitude is observed to decrease during the burst rise while the 
emitting area increases, reaching a constant value during the 
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Table 1 

Properties of Observed Type I Burst Oscillations 
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4U 1636-54 


1996 Dec 28 (22:39:22) 
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RJ 2-A 


Y 


1,2'' 




1996 Dec 29 (23:26:46) 
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ai2 




Si 3 
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1996 Dec 31 (17:36:52)= 
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ss 3 


Y 


2,3'' 
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1997 Jul 26(14:04:19) 


7 


329.85 ±0.1 


2.5 


1.88±0.25 


7.7±0.3 


Y 


4" 




1997 Jul 30(12:11:58) 
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330.55 ±0.02 


1.6 


4.02 ±0.07 


4.8 ±0.3 


N 


4" 


4U 1728-34 


1996 Feb 16 (10:00:45) 
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364.23 ±0.05 


2.4 


3.52 ±0.28 


6.6±0.1 


Y 


4,5" 




1997 Sep 9 (06:42:56) 


7 


364. 10 ±0.05 


2.1 


1.84±0.15 


5.9 ±0.2 


? 
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KS 1731-26 


1996 Jul 14 (04:23:42) 


Y 


524.61!«;» 


0.96 ±0.04 


4 1+1.3 


1.8 


N 


6,7" 




1999 Feb 27 (17:23:01) 
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524.48!0;»^ 


1.06 ±0.05 


2 6+^-^ 
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2.0 


N 


7" 
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1996 May 15 (19:32:10) 
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^ 1 
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8" 


(MXB 1743-29?) 
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1997 Mar 1 (23:27:39) 
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549.76 ±0.04 


00+0.03 


2.69±0.15 


4.3 


N 


9,10" 



"In these cases, the frequency evolution was fitted by an exponential model v{t) = ug-Aue 

''For these bursts, uq is the frequency seen in the tail, and Au how much the frequency changes during the burst. 

'^This burst observed an episode of spin down in the tail (Strohmayer 1999b). 

References.— (1) Strohmayer et al. 1998c; (2) Miller 2000; (3) Strohmayer 1999a; (4) Strohmayer & Markwardt 1999; (5) Strohmayer et al. 1996; (6) Smith, 
Morgan & Bradt 1997; (7) Muno et al. 2000; (8) Strohmayer et al. 1997a; (9) Zhang et al. 1998; (10) Fox et al. 2000. 



decay (Strohmayer, Zhang & Swank 1997b; Strohmayer et al. 
1998c). 

The oscillations during the burst rise are well-explained by 
the picture of a spreading hotspot. However, several myster- 
ies remain. First, oscillations are often seen in the burst tail, 
when the whole surface of the star has presumably ignited. The 
cause of the azimuthal asymmetry at late times is not under- 
stood. Secondly, for those objects with i/o « 550 Hz in Table 
there is evidence that the burst oscillation frequency is twice the 
spin frequency of the star, including observations of a w 290 Hz 
subharmonic (Miller 1999) and the fact that the kHz QPO sep- 
aration is approximately half the burst oscillation frequency in 
these sources (see van der Klis 2000 for a review). What might 
cause an m = 2 azimuthal asymmetry is unknown. Third, oscil- 
lations are not seen in all sources, or all Type 1 bursts from the 
same source. 

An initially puzzling feature of the observations was that the 
oscillation frequency often changes during the burst, increasing 
by A^Ri 1-2 Hz in the burst tail. Strohmayer et al. (1997a)pro- 
posed a simple explanation — that this frequency shift results 
from angular momentum conservation. The slight hydrostatic 
radial expansion (contraction) of the burning layers as the tem- 
perature increases (decreases) results in spin down (spin up) 
if angular momentum is conserved. The time for radial heat 
transport from the burning layers to the photosphere is about 
one second, which means that the layers are puffed up and spun 
down by the time the observer sees the burst. As the burn- 
ing layers cool during the tail of the burst, they contract and 
spin up. Strohmayer & Mai'kwardt (1999) and Miller (1999, 
2000) have modelled the observed frequency evolution. They 
find that, in the evolving frame of the burning shell, the os- 
cillations are coherent, as expected if they are due to rotation. 
Calculations show that the change in thickness of the burning 
layers during a burst is Az « 20 m (Hanawa & Fujimoto 1984; 
AyasU & Joss 1982; Bildsten 1998). A simple estimate of the 
spin down of the burning shell due to this change in thickness is 



then Av ^ v,{2Az/R) « 1 Hz(i^,/300Hz)(Az/20m)(10km/7?) 
where R is the neutron star radius and Vs the spin frequency. 
This roughly agrees with the observed values (Table |l]). 

In this paper, we make a first attempt at understanding the 
evolution of the neutron star atmosphere during a Type 1 X- 
ray burst on a rotating neutron star. We use hydrostatic models 
of the neutron star atmosphere to calculate the expansion and 
resulting spin evolution during a burst. Our aim is to investi- 
gate whether the observed frequency changes during bursts are 
consistent with spin down of the burning shell, thus lending 
support to the interpretation that the burst oscillation frequency 
is intimately related to the neutron star spin. This simple pic- 
ture demands that the hot burning material decouples from the 
cooler underlying ashes and conserves its angular momentum, 
completing a few phase wraps with the underlying star during 
the burst. We thus examine mechanisms that might couple the 
burning layers to the neutron star, and ask whether it is plausi- 
ble that the burning layers can remain decoupled for the « 10 s 
duration of the burst. 

We stress that we consider only one-dimensional models for 
these hydrostatic and coupling calculations. We do not consider 
the complex question of how the burning front spreads around 
the star during the burst rise, what determines the number of 
hotspots on the surface, or what causes the asymmetry at late 
times during the cooling tail of the burst. We leave these ques- 
tions for future investigations. 

We start in §2 by describing the hydrostatic structure of the 
atmosphere during fuel accumulation and the X-ray burst. In 
§3, we calculate the expansion of the atmosphere and show that 
the expected spin down of a decoupled layer is consistent with 
observed values. We consider heat transport through the atmo- 
sphere, and ask how a single coherent frequency is transmitted 
to the observer In §4, we discuss hydrodynamic mechanisms 
that could transport angular momentum within the burning lay- 
ers or couple the hot burning material to the underlying colder 
and denser ashes. In §5, we summarize our results and discuss 
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the many remaining puzzles. Finally, we present our conclu- 
sions in §6. 

2. THERMAL STRUCTURE AND EXPANSION OF THE BURNING 
LAYERS 

Most neutron stars in LMXBs accrete hydrogen and he- 
lium rich material from their companions at rates M ^ 10~" - 
IO'^Mq yr"'. For accretion rates M ^ 2 x 1O"'°M0 yr"' (see 
Bildsten 1998 and references therein), the accumulating hydro- 
gen is thermally stable and burns via the hot CNO cycle of 
Hoyle & Fowler (1965). The temperature of most of the at- 
mosphere is ^ 8 X 10^ K so that the time for a proton capture 
onto a ''^N nucleus is less than the time for the subsequent beta- 
decays. This fixes the energy production rate at the value 



en = 5.8 X 10' 



/ ZcNO \ 

['om 



ergs g s 



(1) 



where Zcno is the mass fraction of CNO nuclei. This energy 
production rate is independent of temperature or density, and 
simply set by the beta-decay timescales of ''^O (half-life 71s) 
and '^O (half-life 122 s). Because the hydrogen burns at a con- 
stant rate, the time to burn all of it depends only on the metallic- 
ity and initial hydrogen abundance. The hydrogen mass fraction 
X in a given fluid element changes at a rate dX /dt = -en /Eh, 
where Eh ~ 6.7MeV/proton « 6.4 x 10"^ergg"' is the energy 
release per gram from burning hydrogen to helium. The time to 
burn all the hydrogen is then 



tn ~ 22 hours 



0.01 



Zcno/ VO-71 



^0 



(2) 



where Xq is the initial hydrogen mass fraction. 

The X-ray burst is triggered when helium burning becomes 
unstable at the base of the accumulated layer, at a density 
^ 10^-10''gcm"-' and temperature w 2 x 10^ K. The composi- 
tion at the base of the layer depends on how much hydrogen has 
burned during the accumulation (Fujimoto, Hanawa & Miyaji 
1981, hereafter FHM; Bildsten 1998), which is determined by 
the local accretion rate, m. The local Eddington accretion rate 
is iTiEdd = ^nipC / {I + X)RaT:h, where crxh is the Thomson scat- 
tering cross-section, nip is the proton mass, c is the speed of 
light, and R is the stellar radius. In this paper, we use the 
Eddington accretion rate for solar composition (Xq = 0.71) and 
R = 10 km, niEAd = 8.8 x 10"* g cm"^ s"', as our basic unit for the 
local accretion rate (for a 10 km neutron star, this corresponds 
to a global rate M = 1.7 x lO'^Mgyr"'). For m ^ 0.03;«Edd, 
the helium burning becomes unstable before all the hydrogen 
is burned and the helium ignites and burns in a hydrogen rich 
environment. At lower accretion rates m <^ 0.03 mEdd, there is 
enough time to burn all the hydrogen, and a pure helium layer 
accumulates which eventually ignites. We now describe sim- 
ple models of the atmosphere immediately prior to (§2.1) and 
during (§2.2) the X-ray burst for both of these cases. At lower 
accretion rates still, lii < O.Ol/WEdd, the hydrogen burning be- 
comes unstable and triggers a mixed hydrogen/helium burning 
flash (FHM). 

2.1. The Accumulating Atmosphere 

The neutron star atmosphere is in hydrostatic balance as 
the accreted hydrogen and helium accumulates. The pressure 
obeys dP/dz = -pg, where p is the density and the gravita- 
tional acceleration g = GM/R^ is constant in the thin enve- 
lope. A useful variable is the column depth y (units g cm"^). 



defined by dy = —pdz, giving P = gy. As the matter accumu- 
lates, a given fluid element moves to greater and greater col- 
umn depth. In this paper, we take gu = g/lQ^* cms"^ = 1.9, 
appropriate for a M = IAMq and R = 10 km neutron star As 
described above, the hydrogen burning rate is a constant so that 
the change of X with column depth y is dX /dy = -en / mEn, 
where we take Eh = 6.4 x lO'^ erg g~' (we have neglected the 
« 0.5 MeV/proton lost as neutrinos in the hot CNO cycle, see 
Wallace & Woosley 1981). Integrating this equation, we find 
the hydrogen abundance as a function of depth is 







(y/yd)] y>yd, 
y<yd. 



(3) 



where Xq is the initial hydrogen abundance, and the column 
depth at which the hydrogen runs out is y^i = XoinEfj /en, or 



yd = 6.S X lO'^gcm" 



O.lmEdd 



/0.01\ ( Xq\ 

\ZcNoJ \o.nJ 



(4) 



If helium ignites at a column depth y < yd, a mixed hydro- 
gen/helium burning flash occurs, otherwise a pure helium layer 
accumulates which eventually ignites at a column y > y^. 

The thermal profile of the accumulating layer is described by 
the heat equation, 

dT_ _ 3kF 
dy AacT^ ' 

where k is the opacity and F is the outward heat flux. The en- 
tropy equation is 



dF dT cpTrh / dlnT 
^ dy ''^ dt y \ dlny 



(6) 



(Bildsten & Brown 1997), where n = (dlnT /d\nP)s, cp is the 
heat capacity at constant pressure, and the terms on the right 
hand side describe the compression of the accumulating matter. 
The compressional terms contribute « cpT « ShgT /Ijimp « 
20 Ts keV per accreted nucleon, where p, is the mean molecular 
weight, and = T/IO** K. In comparison, the hot CNO cy- 
cle hydrogen burning gives « IXq MeV per accreted nucleon 
for pure helium ignition, or « 7Xo(yh/yd) MeV for mixed H/He 
ignition, where yb is the column depth at the base. There is 
additional flux from heat released by electron captures and py- 
cnonuclear reactions deep in the crust, giving « 100 keV per 
nucleon (Brown & Bildsten 1998; Brown 2000). For the pur- 
poses of this paper, we neglect the compressional terms in the 
entropy equation, and take 



(7) 



but include a flux at the base (y = yb) of Fh= 150 keV per nu- 
cleon as an approximation to the heat from the crust and com- 
pressional heating. 

We find the temperature profile by integrating equations (||) 
and (^ from the top of the atmosphere to the base at col- 
umn depth yh- At the top, which we arbitrarily place at y = 
5 X 10"^gcm"^, we set the temperature using the analytic radia- 
tive zero solution for a constant flux atmosphere with Thomson 
scattering opacity (the solutions are not sensitive to this upper 
boundaiy condition). For y < we take e = en, and for y >yd 
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Table 2 

Plane Parallel Ignition Conditions* 





ZcNO 


r 


y 


X 


y 


P 


yjirt 


Az(90%)'' 






(10" K) 


(lO^'gcm-^) 






(lO'-gcm-^) 


(h) 


(m) 


Pure He Ignition 


0.01 


0.005 


1.49 


5.40 


0.0 


0.995 


2.28 


170 


4.2 




0.01 


1.35 


12.9 


0.0 


0.99 


4.14 


407 


5.1 




0.02 


1.22 


34.0 


0.0 


0.98 


8.01 


1073 


6.9 


0.015 


0.01 


1.56 


3.67 


0.0 


0.99 


1.75 


11 


3.8 




0.02 


1.41 


8.50 


0.0 


0.98 


3.11 


179 


4.4 


0.02 


0.01 


1.74 


1.92 


0.0 


0.99 


1.12 


30 


3.8 




0.02 


1.57 


3.62 


0.0 


0.98 


1.74 


57 


3.5 


0.03 


0.02 


1.83 


1.56 


0.0 


0.98 


0.97 


16 


3.6 


Mixed H/He Ignition 


0.015 


0.005 


1.73 


2.05 


0.01 


0.99 


1.16 


43 


4.2 


0.02 


0.005 


1.78 


2.18 


0.15 


0.85 


1.04 


34 


4.5 


0.03 


0.005 


1.86 


2.35 


0.31 


0.69 


0.94 


25 


4.9 




0.01 


1.89 


1.67 


0.14 


0.85 


0.87 


18 


4.2 


0.1 


0.005 


2.05 


2.70 


0.57 


0.42 


0.84 


8.5 


5.6 




0.01 


2.13 


2.04 


0.50 


0.49 


0.71 


6.4 


5.1 




0.02 


2.20 


1.51 


0.40 


0.58 


0.62 


4.8 


4.7 


0.3 


0.005 


2.24 


2.62 


0.67 


0.33 


0.76 


2.8 


5.8 




0.01 


2.33 


2.14 


0.64 


0.35 


0.66 


2.3 


5.5 




0.02 


2.44 


1.70 


0.59 


0.39 


0.57 


1.8 


5.2 



"Conditions at the base of the accumulated column at ignition. We take M = 1.4M0, R = 10km, 
and g = GM/R =1.9x10 cms . The ignition pressure is P = = 1.9^8 x 10 ergcm and the 
accumulated mass over the whole surface is AM = AirR^y = 1.3y8 x 10^' g. We take the flux from 
deeper in the star to be F4 = 150 keV per nucleon. 

''mEdd = 8.8 x lO^gcm^^s"', equivalent to a global rateMsdd = 1.7 x lO^^Moyr"'. 

■^The time to accumulate the unstable column. 



""The height above the base which contains 90% of the mass. 
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we set e = 0. The flux at the top F, is set by the energy re- 
lease from hot CNO burning in the atmosphere and the flux at 
the base, giving F, = Fb + enyh or Ft = Fb + enyd, whichever is 
smaller^. The opacity k has contributions from electron scat- 
tering, free-free absorption and electron conduction, which we 
calculate as described in Schatz et al. (1999). Since in the hot 
CNO cycle the seed nuclei are mostly '"^O or '^O waiting to f3- 
decay, we take the gas to be a mixture of hydrogen (mass frac- 
tion X given by eq.[|^]), '"'O and '^O (mass fraction Zcno) and 
helium (mass fraction Y = 1 -X -Zcno)- The ratio by number 
of ''^O to '^O is given by the ratio of the beta-decay timescales, 
giving Zi4 = 0.352ZcNO and Z15 = 0.648Zcno- To obtain a 
simple analytic estimate of the base temperature, we first in- 
tegrate equation (^ with en a constant to find F(y) = eniyh-y), 
where we ignore compressional heating and the flux from the 
base. Substituting this into equation (||) and integrating assum- 
ing constant opacity, we find T'^iy) = {2>hi/ ac)eHyyb{^ -y/'^yt), 
giving 



Tb « 2.6 X 



lO'^K 



0.08 cm2 g- 



1/4 



yb 



Zc 



NO 



0.01 



2 X lO'^gcm" 



1/4 



1/2 



(8) 



where we have inserted a typical value for k. 

We find the column depth at the base of the accumulated 
layer just before the thermally unstable helium ignition by com- 
paring the temperature sensitivity of the heating and cooling 
rates (FHM; Fushiki & Lamb 1987b; Bildsten 1998). The heat- 
ing rate due to the triple-alpha reaction is 



and composition at the base, the time to accumulate the unsta- 
ble column yb jm, and the physical distance from the base to the 
place where y = yb/lO, Az(90%). For mixed H/He (pure He) 
ignition, Az(90%) w 5 m (3.5-7 m). 

Figure |l](a) shows the temperature profile of models with 
mixed H/He ignition. The solid lines are for Zcno = 0.01 and 
(bottom to top) m/m-EdA = 0.03,0.1,0.3. The dashed lines are 
for m = 0.1 mEdd and Zcno = 0.005 (bottom curve) or 0.02 (top 
curve). Figure ^b) shows temperature profiles for pure He 
ignition. The solid lines show Zcno = 0.01 and (bottom to 
top) m/niEdd = 0.01,0.015, and 0.02. The dashed lines are for 
m = 0.01 mEdd and Zcno = 0.005 (top curve) or 0.02 (bottom 
curve). The black dots show the depth where hydrogen runs 
out (at a column y^, eq. |Q|). 

For the mixed H/He ignition models, the ignition temperature 
and density do not depend sensitively on m. Slight differences 
arise because at higher accretion rates less helium is made from 
hydrogen burning, requiring a higher density and temperature 
for ignition. For pure He ignition, the ignition column and tem- 
perature are much more sensitive to accretion rate. After the 
hydrogen is burned, there is a slight temperature gradient to 
carry the flux from the base Fb, but the temperature at igni- 
tion is mainly set by the temperature at the base of the hydro- 
gen burning shell (FHM; Wallace, Woosley and Weaver 1982). 
This temperature is greater at higher accretion rates, leading to 
a smaller column depth at ignition. 

The effect of increasing metallicity for mixed H/He ignition 
models is to increase the energy generation rate, and thus the 
temperature at the base. Together with the increased rate of 
helium production, this allows He ignition at a lower column 
depth, as seen in Figure [ij a). For pure He ignition, the effect of 
increasing metallicity is exactly opposite. The ignition temper- 
ature is lowered with increasing metallicity because the hydro- 
gen runs out at a smaller column, giving a smaller temperature 
at the base of the hydrogen burning shell. Ignition then requires 
a greater column depth at the base of the almost isothermal pure 
helium layer, as seen in Figure |l](b). 

The m at which the transition between pure He and mixed 
H/He ignitions occurs depends on metallicity. The m = 
0.02 mEdd case in Figure [] (right panel) just burns the hy- 
drogen before ignition; this is the transition m for a metallic- 
ity Zcno = 0.01. For Zcno = 0.005, the transition occurs for 
m « 0.01 njEdd, and m > 0.03 mEdd for Zcno = 0.02. This 

13/18 

is in rough agreement with the oc Z^^j^'q scaling estimated by 
Bildsten (1998). The transition accretion rate increases with in- 
creasing metallicity because the time to burn all the hydrogen 
decreases, allowing a pure helium layer to build up even at large 
accretion rates. 

Our results agree well with those of previous authors. The 
calculations we have presented here are similar to those of Taam 
(1980) and Hanawa & Fujimoto (1982). Hanawa & Sugimoto 
(1982) simulated pure He ignition bursts on neutron stars with 
g\4 = 0.93 and 7.1. Conditions at the base of the hydrogen burn- 
ing shell that they find agree with our calculations to 10%, and 
the helium ignition column agrees to within a factor of two. We 
find similar agreement with the pure He ignition calculations 
of Wallace, Woosley & Weaver (1982, hereafter W82). How- 
ever, the density at the base of the hydrogen burning shell that 

^As pointed out by Taam & Picklum (1978), some helium burning occurs during accumulation, which increases the number of CNO nuclei at ignition (for exam- 
ple, see Fujimoto et al. 1987, Figure 4). FHM and Hanawa & Fujimoto (1982) estimated the increase in CNO abundance, concluding that the effect on the ignition 
temperature and density would be small because of the large temperature sensitivity of the triple alpha reaction. In calculations which include helium burning during 
accumulation, we find that helium burning reactions decrease the ignition column by 10-20%, increase the ignition temperature by less than a few percent and increase 
the metallicity at ignition by factors of 2-3. For the purposes of this paper, we adopt simple models with hot CNO burning only. 



£30 



= 5.3 X 



lO^'ergg- 



2v3 



exp 



-44 



(9) 



where / is the screening factor (Fushiki & Lamb 1987a). In ad- 
dition, '^C rapidly captures protons once it is made, increasing 
the energy release from the triple alpha reaction. To account 
for this, we multiply £30, by a factor 1 +2i2/G3a = 1-9, where 
Q3a = 7.274 MeV, Qn = QC^C + 2p -^'^ q) = 6.57 MeV. An 
effective local cooling rate is obtained from equations (|]) and 
(0), giving 

iKy^ 

As the matter accumulates, the column depth at the base yb in- 
creases until deia/dT > decooi/dT, at which point a thermal 
instability occurs (FHM). We find the temperature profile at 
ignition by choosing yb such that dey,a/dT = de^oox/dT at the 
base. In Figure |l], the hatched region in the temperature-column 
depth plane indicates where helium ignition occurs for abun- 
dances ranging from the initial value (Y = 0.3) to pure helium 
(F = 1.0). 

The conditions at the base of the accumulated column at the 
time of ignition for a variety of metallicities and accretion rates 
are shown in Table ^ We have separated our solutions into 
two groups depending on whether the helium unstably ignites 
before (mixed H/He ignition) or after (pure He ignition) the hy- 
drogen is completely burned. We give the temperature, density 
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Table 3 

Models with a Convection Zone 



{yb-yc)lyb'' 




Azc Az(90%) 


Pb 




X 


Y 






(10' K) 


(m) (m) (10 


' g cm-') 












Mixed Ignition 














{m = Q.i riiEAi; ZcNO 


= 0.01; yfc = 2.04 X 10** g 


cm ; 


Az(yU7c ) = 


5.1m) 






0.99 


1.7 0.89 


61.5 40.9 


0.83 


0.46 


0.60 


0.39 


0.67 


0.95 


1.7 1.15 


48.5 40.7 


0.84 


0.46 


0.60 


0.39 


0.67 


0.8 


1.5 0.74 


18.3 24.4 


1.41 


0.67 


0.58 


0.41 


0.68 


0.5 


1.2 0.31 


5.5 15.9 


2.30 


0.86 


0.55 


0.44 


0.70 






Pure He Ignition 














{lii = 0.015 OTEdd; ZcNo = 0.01; y,, = 3.67 x 10** 


gcm- 


; Az(90%) -- 


= 3.8m) 






0.99 


2.0 1.06 


45.6 30.2 


2.04 


0.42 


0.09 


0.90 


1.16 


0.95 


1.95 1.21 


29.8 25.0 


2.45 


0.45 


0.07 


0.92 


1.20 


0.8 


1.6 0.59 


8.9 13.1 


5.16 


0.76 


0.01 


0.98 


1.32 


0.5 


1.3 0.35 


3.1 10.0 


7.38 


0.90 


0.00 


0.99 


1.35 



"The fraction of the accumulated mass that becomes convective. 
''The ratio of gas pressure to total pressure at the base. 
'FEdi = 1.67 X 10-^ erg cm"^ s"'. 



we compute is twice as great as that given in Table 1 of W82, 
whereas the pressure and temperature agree. We do not know 
the reason for this discrepancy, but it may be that W82 used a 
fixed solar composition for the hydrogen burning shell, giving 
a different value for /i^ at the base. Hanawa & Fujimoto (1984) 
simulated mixed H/He ignition bursts. Taking into account their 
gravity, gi4 = 3.3, our ignition column and temperature agree 
with their calculations to 10-30%. 

2.2. The Atmosphere During the Burst 

Many authors have simulated X-ray bursts for both mixed 
H/He and pure He ignitions, for reviews see Lewin, van 
Paradijs, & Taam (1993), and Bildsten (1998). In this section, 
we make simplified models which allow us to calculate the hy- 
drostatic expansion of the atmosphere during a burst. We do not 
discuss the lateral spreading of the burning front in this paper 
In addition, our models are not appropriate for the radius ex- 
pansion phase of bursts with super-Eddington luminosities (ra- 
dius expansion bursts; Lewin, van Paradijs, & Taam 1993). We 
discuss the effe ct o f radius expansion on the observability of 
oscillations in §3.5 We start in §2.2.1 by considering the con- 
vective stage of a burst. As we discuss later, convection is im- 
portant because it likely enforces rigid rotation, and m ay affect 
the observability of a coherent signal (§^. In §2.2.2 



pute fully-radiative atmospheric models, appropriate for those 
bursts that do not become convective, or for the later stages of 
a convective burst, when the convection zone retreats. 

2.2.1. Models with Convection 

In one-dimensional models, the energy release from the tem- 
perature sensitive helium burning reaction makes the early 
stages of many bursts convective. The convection zone expands 
outwards, extending over a few scale heights to the radiative 



outermost layers (Joss 1977). For example, Hanawa & Fuji- 
moto (1984) computed a mixed H/He ignition model, in which 
the atmosphere was convective for wis after ignition, reaching 
a maximum extent of « 99% of the accumulated mass. Hanawa 
& Sugimoto (1982) found for pure He ignition that a convec- 
tion zone rapidly grew to encompass most of the atmosphere, 
and then shrunk back, disappearing « 1 s later when the nuclear 
fuel was almost exhausted. 

The convection in the neutron star atmosphere is very effi- 
cient, since the time for sound to cross a scale height {cjg ^ 
10~''s) is much shorter than the local thermal time (~ 1-lOs). 
Thus the atmosphere has a nearly adiabatic profile when con- 
vective. For a given temperature Ti, and column depth yi, at the 
base, the thermal profile of the convective zone just follows the 
adiabat T = Tb(y/yb)", where « = {dlnT /dlnP)s. We take yh to 
be the value given by the settling solution at ignition; the con- 
vection zone profile is then determined by the single parameter 
Th when the mean molecular weight, /i, is fixed. 

During the burst, the temperature at the base reaches ;> 10^ K, 
high enough that radiation pressure is important (Joss 1977). At 
these temperatures, the degeneracy is lifted (we find the elec- 
tron pressure differs from the ideal gas value by less than a few 
percent) so the total pressure at the base is 



Pb- 



PbkBTj, 1 A 
^inip 3 



(11) 



where pb is the density at the base. However, the total pres- 
sure is also fixed by the weight of the overlying atmosphere, 
which gives Pt = gyt- Thus Tf, cannot exceed the critical value 

Tlnax = Ogyh/af^'', or 



r„,, = 2.2x lO-'K 



3 X 10**gcm"2 



1/4 



1/4 



(12) 
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Table 4 
Fully Radiative Models 



F 1 Fedd 7* 


Az(90%) 


Pb 


/3;, Notes 


(10' K) 


(m) 


(10^ gem" 








Mixed Ignition 






(m = 0. 1 mEdd ; ZcNO 


= 0.01;>'6 = 


2.04 X 10** I 




J. 1 m) 


1.15 1.49 


29.7 


1.53 


0.68 


a 


1.29 


25.7 


2.13 


0.82 


b 


0.74 1.35 


22.3 


1.94 


0.78 


a 


1.17 


20.1 


2.50 


0.88 


b 


0.31 1.12 


15.6 


2.67 


0.90 


a 


0.96 


14.5 


3.20 


0.94 


b 




Pure He Ignition 






(m = 0.015 OTEdd; ZcNO = 0.01; jj, 


= 3.67 X 10** 


gcm^^; Az(90%) 


= 3.8m) 


1.21 1.57 


16.9 


5.44 


0.78 


a 


1.38 


15.4 


6.83 


0.87 


b 


0.59 1.34 


11.9 


7.10 


0.88 


a 


1.17 


11. 1 


8.35 


0.93 


b 


0.35 1.19 


9.9 


8.22 


0.93 


a 


1.04 


9.4 


9.39 


0.96 


b 



"Constant flux. 
''Constant e. 



As Th approaches Tmax, radiation pressure becomes increasingly 
important, forcing the density to decrease. We show below that 
this greatly enhances the expansion of the atmosphere. 

As clearly argued by Joss (1977), the convection zone can- 
not extend all the way to the photosphere. A radiative layer is 
always needed to transport the heat to the photosphere. We as- 
sume that most of the burning occurs in the convective layer, in 
which case the radiative atmosphere carries a constant flux and 
is described by the heat equation (^ with constant F . For the 
convective zone, we integrate dT /dy = n{T /y) for a given base 
temperature Ti, (the super-adiabaticity needed to carry the flux 
in the convective zone is negligible when the convection is so 
efficient; see Cox & Giuli 1968). For these simple models, we 
treat the flux F and base temperature Th as free parameters (in 
reality, the temperature throughout the convective region deter- 
mines the nuclear energy generation rate and thus the flux). To 
compute the thermal profile of the atmosphere, we integrate the 
convection zone outwards from the base, and the radiative zone 
inwards from the surface, varying the column depth of the top 
of the convection zone yc until the temperature matches at the 
interface. 

We compute convective models for a range of fluxes and base 
temperatures for two of the settling solutions of Figure [l|. The 
first is for the mixed H/He ignition model with m = 0. 1 mEdd and 
ZcNO = 0.01. Figure ^a) shows temperature profiles of con- 
vective models with flux F/Fem = 0.31,0.74,0.89 and 1.15, 
the fraction of the accumulated matter that becomes part of 
the convective zone, {yh-yc)/yh = 0.5,0.8,0.99 and 0.95, and 
base temperature ri/10^K= 1.2,1.5,1.7 and 1.7 respectively. 
Here we measure flux in units of the Eddington flux (for the 
accreted solar composition) Z^Edd = GMniEdd/R, giving /^Edd = 



1.67 X 10^^ergcm"^s"'. Figure ^b) shows convective mod- 
els coiTesponding to the pure He ignition settling solution with 
m = 0.015 wiEdd and Zcno = 0.01. We make convective mod- 
els with flux F/FEdd = 0.35,0.59,1.06 and 1.21, the fraction 
of the accumulated matter that becomes part of the convective 
zone, (yb-yc)/yb = 0.5,0.8,0.99 and 0.95, and base temperature 
ri/10'*K= 1.3,1.6,2.0 and 1.95 respectively. Properties of aU 
these models are given in Table |[ In each case, we take the 
composition profile of the atmosphere to be that immediately 
prior to ignition, but assume full mixing of matter in the con- 
vective zone (we average over the convection zone according to 
Xi = J dyXi{y)/{yh-yc) for each species /). 

We give the thickness of the convection zone Azc in Table ^ 
as well as the height above the base which contains 90% of the 
mass, Az(90%). A simple analytic estimate of the thickness, 
Azc, is obtained as follows. First, the adiabatic index of the gas 
is 

^{d}nT\ 8-6/3 
"-\d\np)^ 32-24/3-3/32' 

(Clayton 1983), where the ratio of gas pressure to total pressure 
is 

P=1-Q,QU(^X( X( — V- (14) 

^ yiO'^Kj V3xl08gcm-2y Vl.9/ 

For /3 = 1 (P = Pg), we obtain the ideal gas value n = 2/5; for 
/3 = (P = P,), n= 1 /4. We take P oc T", and presume /3 is a 
constant throughout the convective zone, so that n is a constant, 
given by equation (|T3|). Integrating dz = —dP/ pg, we find the 
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thickness of the convection zone is 



3. SPIN EVOLUTION OF THE BURNING LAYERS 



n 



l-(^ 
Ph 



(15) 



when it extends to a pressure P,. We denote the scale height at 
the base as Hb = Pb/ Pbg = kBTb/fJ.mpgf3 and assume ^ is con- 
stant, giving 



fJ-P \ n J \gi4 



l-(^ 
Pb 



(16) 



Using the value of P at the base (denoted (3b in Table to evalu- 
ate equation (|6|), we find this analytic estimate agrees to better 
than 5% with the numerical integrations. Equation (|6|) shows 
that radiation pressure strongly affects the radial extent of the 
convective zone, Az^ oc 1 /f3. 

We show the density as a function of height above the base 
for the mixed H/He and pure He ignition models in Figure |^(a) 
and (b) respectively. The jump in density at the interface be- 
tween the convective and radiative zones is due to the jump in 
composition because the convective zone is fully mixed. The 
black dots mark the height which encloses 90% of the mass as 
found numerically. This ranges from w 15-50 m for the mixed 
H/He models, and « 10-30 m for the pure He models (Table 
compared to w 5 m before the burst (Table The difference 
between the mixed H/He and pure He models arises because 
Azc depends strongly on the mean molecular weight (eq. JT^ ] 
gives Azc oc l//i). For a mixture of H and He, /i = 4/(8X+3}'), 
giving 11 = 4/3 for pure He, but /i = 0.6 for solar composition, 
almost a factor of two different. The mean molecular weight 
for each model is given in Table |[ 

2.2.2. Fully-Radiative Models 

We consider two different kinds of radiative models: a radia- 
tive atmosphere carrying constant flux, and a radiative atmo- 
sphere with a constant energy production rate e. The first case 
is relevant when the energy production is concentrated near the 
base, for example in the initial He burning stages of a burst 
which does not convect. The second case is relevant when the 
burning region is more extended. For example, in bursts with 
mixed H/He ignition, the helium burns rapidly (by purely strong 
reactions) during the initial stages of the burst, whereas the hy- 
drogen burns more slowly (by the rp-process involving weak 
reactions; Wallace & Woosley 1981), giving rise to the long 
tails seen in some X-ray bursts (for a recent example, see Fig- 
ure 2 of Kong et al. 2000). Because the timescale for H burning 
by the rp process (^ 10-100 s) is long compared to the time for 
radiative heat transport {^^ 1-10 s), the atmosphere is radiative 
during the rp process tail of the burst (Hanawa & Sugimoto 
1982). 

We give details of our fully-radiative models in Table ^ For 
the constant e models, we choose e to give the required flux 
at the top, F = eyb+Fh- We assume the composition of the at- 
mosphere is the same as immediately prior to ignition. Figure 
Q shows the temperature profiles. We compute models with 
F/Fem = 0.31,0.74 and 1.15 (mixed ignition) and F/FEdd = 
0.35, 0.59 and 1 .21 (pure He ignition). Figure || shows the den- 
sity profiles for these models. For the mixed ignition models, 
Az(90%) « 15-33 m; for the pure He models, Az(90%) w 9- 
17 m. Again, the pure He models are less extended because of 
the larger fi. 



We have shown that the atmosphere expands outwards by 
« 10-50 m (« 5-25 m) during a mixed H/He (pure He) burst. In 
this section, we first calculate how the spin of the burning layers 
evolves during the burst, assuming they conserve angular mo- 
mentum as they hydrostatic ally expand and contract, and that 
they remain decoupled from the bulk of the neutron star (§3.1- 
3.3). In §3.4, we point out (in agreement with Miller 2000) that 
for the observer to see a coherent oscillation requires either a 
mechanism to enforce rigid rotation of the burning layers (such 
as convection), or that the burning layers be geometrically thin. 
In addition, we consider the heat propagation through the radia- 
tive layers to the surface, and discuss how this can "wash-out" 
any coherent pulse emanating from deeper layers. We conclude 
in §3.5 by discussing what happens during radius expansion 
bursts. 

3.1. Expansion and Spin Down 

We assume that the action of hydrodynamic instabilites (Fu- 
jimoto 1988, 1993) or a weak poloidal magnetic field will force 
the accumulating pile of fuel to be rigidly rotating with the 
spin frequency of the star, The time for angular momen- 
tum to diffuse across a scale height H due to molecular vis- 
cosity I' alone is t^.jse ~ /v. In §5, we show that this time is 
tvisc « 38h (i//5 m)^ for Tg « 2, comparable to the time between 
bursts. Thus, even microphysical mechanisms might bring the 
accumulated fuel into corotation. 

We further assume that the burning is spherically symmetric, 
ignoring the complicated, and important, question of how the 
burning spreads over the stellar surface, and the fact that some 
asymmetry is needed to give an observable oscillation. Thus we 
assume the angular velocity is constant on spherical shells, and 
depends on radius only. This remains true as the atmosphere 
expands, since a rigidly rotating spherical shell stays rigidly ro- 
tating as it expands outwards (because the fractional change in 
the distance from the rotation axis is the same for all latitudes). 
Conservation of angular momentum demands that p-Vl{r) re- 
mains constant as a fluid element expands outwards, giving the 
spin frequency as a function of column depth during the burst 



^ \R+zi{y)) ' 



(17) 



where zi {y) {ziiy)) is the height of column y above the base be- 
fore (during) the burst. We take the base at yb to be at fixed 
radius R, even though the old ashes below heat up and expand 
a little. We ignore the general relativistic corrections to the fre- 
quency shift, which Strohmayer ( 1 999a) estimates are a « 10%- 
20% effect. 

Equation ( |l7| ) shows that the spin frequency depends on 
depth in the radiative layers. The picture is different in a con- 
vective zone. We expect the convective motions will rapidly 
bring the convection zone into rigid rotation, since the convec- 
tive turnover time is very short (^ 10"^ s). Thus, if there is no 
angular momentum transport between the convective zone and 
neighboring layers, the spin frequency of the convective zone is 



(18) 



where /] is the total moment of inertia of the material between 
yc and yb before the ignition, and 7^ is the moment of inertia 
of the convection zone. We calculate the moment of inertia 



9 



as follows. A spherical shell at radius r, mass dm = 4nr^pdr 
has a moment of inertia dl = (87r/3)p(r)r^t/r, where the Stt/S 
comes from integrating over angles. Since pdr = -dy, this is 
dl = -{%TT /'i)r'^dy. Thus the moment of inertia is 



i{y) = 



dy[l + 



R 



(19) 



We integrate this equation from the starting values of flc and 
X, towards x ^ 1 as the convection zone vanishes. There is one 
difficulty, however, which is that Az^ and 7^ are functions of not 
only X but also the scale height at the base Hi,. If we assume the 
composition p does not change, then Hi, only depends on Ti,. 
We model Thix) as follows. 



where z(y) is the height above the base, given by integrating 
dz/dy = -\ / p. We take P = gy during this integration, neglect- 
ing the small variation of g with depth, which has an 0(H/R)^ 
effect on the moment of inertia. 

We derive a simple analytic estimate of the moment of iner- 
tia of the convective zone, by assuming /3 = Pg/P is constant, 
so that is given by an equation similar to equation (|l5l). 



Tb(x) = TQ + AT 



l-Xi 



(22) 



Inserting this into equation writing (1 +z/R)'^ 
and integrating, we find 



. = -AM,R^ 
3 



1 + 



4///, 
nR 



1 l-x"- 



1 +4z/R 



(20) 



where x = yc/yb, and AMc = 4TrR^(yt,-yc) is the mass of the 
convection zone. The moment of inertia has two pieces. The 
zeroth order piece, 2AMcR^/3 is the moment of inertia of mass 
AMc concentrated at radius R. The second term is the 0(H/R) 
correction to this because the envelope is extended. Equation 
( ^0| ) reproduces the numerical results to a few percent. 

The lower panels of Figure ^ show A{l/il{y) = (il* - 
^l(y))/il^, as a function of column depth for the different con- 
vective models in Table 1} We define Afl/n so that a positive 
value means that the layers have spun down. Figure Q shows 
Ail/il{y) for the fully radiative models of Table Q The val- 
ues of Ail /n that we find are similar to the observed frequency 
shifts during bursts (Table 1). As we described above, in radia- 
tive zones the spin down depends on height, whereas convective 
zones are rigidly rotating. In the convective models, this gives 
rise to jumps in Aft/fl at the boundaries between the convec- 
tion zone and both the underlying ashes and the overlying ra- 
diative layer. For both convective and fully-radiative models, 
the An/il(y) is smaller by roughly a factor of two for pure He 
ignition as opposed to mixed H/He ignition. This is because the 
pure He models expand outwards less, as we discussed in §2.2. 

3.2. Shrinking of the Convection Zone 



where Tq is the final temperature after the convection vanishes 
(x ^ 1), T() + AT is the initial temperature at the base of the 
convection zone, and x,- is the initial value of x. In reality, x is 
determined by the instantaneous values of F and Tb, so if we 
knew F(Tb) we could map this onto Tbix). Instead, we spec- 
ify Tb{x) then as a check that our choice is reasonable, we can 
compute FiJb). 

Figure ^ shows the results for the model with lii = O.bnEdd, 
•ZcNO = 0.01 (mixed H/He ignition), with initial extent of 
the convective zone (yh-yc)/yb = 0.8, and base temperature 
r/,/10^K= 1.5. As the convection zone retreats, we assume 
the base temperature falls to r/,/10^K= 1.0. We show Afl/Cl 
profiles for the initial model (80% of mass convective), an inter- 
mediate case (50%) and when the convective zone has almost 
vanished (5%). These models have r^/lO^K = 1.5, 1.31, and 
1.03 and F/FEdd = 0.74, 0.49, and 0.21 respectively. The jumps 
in Ail/il at the top and base of the convection zone persist as 
it shrinks. In addition, the matter that becomes radiative spins 
up as the physical thickness of the underlying convection zone 
decreases. This gives rise to an inversion in Aft/fl, as shown in 
Figure ^. One might worry that such an inverted profile would 
be unstable to axisymmetric perturbations. However, all the 
profiles in Figure n have increasing specific angular momen- 
tum with radius. I n ad dition, the Brunt-Vaisala frequency is so 
large, A^^ ^ il? (§4.1 ), that we expect the Solberg-Hoiland cri- 
terion for stability (see e.g. Endal & Sofia 1978) will always be 
satisfied. 



As we described in §2.2.1 



3.3. Spin Up in the Cooling Tail 

As the atmosphere cools in the tail of a burst, its thickness 
decreases, and it spins back up. We compute simple models of 
the cooling atmosphere by assuming it carries a constant flux 
and has a fixed uniform composition. In Figure ^ we show the 
thickness Az(90%) as a function of the flux F /Fem- The solid 
curve is for the m = 0.1 mEdd, Zcno = 0.01 mixed H/He ignition 
model, for which we assume a composition of ^^Kr, since the 
rp process makes elements beyond the iron group (Schatz et al. 
1998 and references therein; Koike et al. 1999). The dashed 
lines are for pure He ignitions, with Zcno = 0.01, m = 0.01 mEdd 
(upper curve) and m = 0.015 ;7jEdd (lower curve), and for a com- 
position ^^Ni. In each case, we show Az(90%) just before igni- 
tion by a horizontal solid line. 

Figure shows that in the limit F/^Edd 0, the thickness of 
the ashes is about half the pre-burst thickness for mixed H/He 
ignition, but similar to the pre-burst thickness for pure He ig- 
nition. This is because electrons (which provide most of the 
pressure support) are consumed in hydrogen burning, but not 

^Since Ic depends on both Ht and x, one can write dlc/dx = die / dx + (dHi, / dx)dlc / dHt- In fact, the first term in equation (Ell) exactly cancels the dic/dx piece 
of the second term, giving dSlcjdx = (Qc /!c)(.dHi,/dx)dIc /dHi,. Thus if the convection zone shrank at fixed Hi,, it would not change its spin frequency because the 
angular momentum taken away by the decoupled fluid exactly cancels the change in the moment of inertia due to the change in x. 



the atmosphere does not remain 
convective for the whole duration of the burst. Eventually, the 
convection zone shrinks back. What is the profile of Aft /ft left 
by the retreating convection zone? Imagine that the convection 
zone retreats by an amount dy^- (so the top of the convection 
zone moves from column to yc+dyc). The material which de- 
couples from the convection zone takes away an amount of an- 
gular momentum given by dJc = -^Tr(R + Azc)*^cdyc/'i, where 
i7c is the spin frequency of the convection zone, and Azc is 
the height of the top of the convection zone above the base. 
The loss of angular momentum from the convection zone must 
be compensated by a change in its spin or moment of inertia 
dJc = d{IcSlc) = Icdftc + ^cdlc- Dividing by dye, and changing 
variables to x = yc/yb, we obtain the differential equation-' 



dQc 
dx 



R 



1 + - 



dh 
dx 



(21) 
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during helium burning. In the tail of a burst, when, for example 
F /Feaa ~ 0.1, the thickness of the atmosphere is w 1 m dif- 
ferent from the pre-burst thickness in both cases. Thus if the 
burning layers remain decoupled as they cool, we expect the 
spin frequency in the tail to be different from the stellar spin by 
a part in « 10 km/1 m « lO'*. We discuss the implications of 



this result in §5.3 



3.4. Heat Transport and Coherence of the Oscillation 

We now turn to the vertical transport of heat, and its effect 
on the coherence and amplitude of the burst oscillations. Fig- 
ures ^ and Q show that in a radiative zone the magnitude of 
the spin down depends on depth. However, a single coher- 
ent frequency is observed during bursts. This implies that ei- 
ther convection or some other mechanism (see §^ is enforc- 
ing rigid rotation in the burning region, or, as pointed out by 
Miller (2000), the burning layers must be geometrically thin. 
Since the fractional change in spin across a height Az in a ra- 
diative zone is w 2Az/R, an observed coherence Q demands 
a burning layer thickness R/2Q. Strohmayer & Markwardt 
(1999) found Q « 4000 for bursts from 4U 1728-34 and 4U 
1702-43 when the observed frequency evolution was accounted 
for, while Muno et al. (2000) found Q « 5000 for burst oscil- 
lations from KS 1731-26. These results imply a burning region 
thickness <^ 100 cm. This is perhaps more likely for pure He 
bursts, because of the temperature sensitivity of He burning re- 
actions. 

Also important is heat transport across the differentially- 
rotating atmosphere. The time for radial heat transport from 
the burning layers to the photosphere can be estimated from the 
entropy equation, cpdT/dt = dF / dy, together with equation (p|) 
for the flux F, which gives a timescale fthem = 3KCpy^/4-acT^. 
In our numerical calculations, we calculate the heat capacity at 
constant pressure cp exactly. We make an analytic estimate for 
a mixture of ideal gas and radiation using cp = {5kB/2iJ,mp)f(f3) 
where /(/3) = {32-24f3-3P^)/5(3^ (Clayton 1983), giving 



ftherm = 0.6 S 



m 

^1 V0-08cm2g-i 



(23) 



The factor /(/?) is unity for (3=1 (P = Pg) and grows in- 
creasingly larger as /? approaches zero (P = Pr), for example 
/(0.85) = 2.6 whereas /(0.4) = 27 (this is because the internal 
energy of a photon gas at constant pressure is independent of 
temperature). Thus at a fixed pressure, the thermal time first 
decreases with increasing temperature (since ftherm oc 1 /T^) but 
starts increasing again once radiation pressure becomes impor- 
tant. 

Equation ( ^3| ) shows that the thermal time (~ 1 s) is much 
longer than the time for hydrostatic readjustment (^ 10"'' s), ex- 
plaining why spin down is not seen in the beginning of a burst. 
By the time the heat released by nuclear burning reaches the 
observer, the layers have expanded and spun down. 

The burning layers revolve around the underlying star in a 
time Ai^~^ = 27r/Af2 wis, comparable to the thermal time. 
The ratio of the shearing time to the thermal time is important. 
If the heat diffuses quickly, ftherm <C 1, the atmosphere sim- 
ply transmits the burning pattern from below, much like pass- 
ing a flashlight behind a piece of paper. If, on the other hand, 
heat diffuses slowly, the oscillations will be smeared out. In the 
limit ftherm Ai^ ^ 1, the atmosphere is heated uniformly by the 
quickly revolving burning layers. 



Figure |] shows ftherm (solid line) and Aj/~' (hatched region 
for = 300-600 Hz) as a function of depth for the convec- 
tive models of Table ^. For mixed H/He and pure He ignitions, 
the wrap around time Az/"' increases for decreasing base tem- 
perature and physical thickness. The thermal time, however, 
decreases with decreasing temperature at a given depth because 
radiation pressure (eq. [|3||) is less important. Figure || shows 
that the mixed H/He ignition bursts have ftherm A f^ « 1, whereas 
the pure He ignition bursts have ftherm Af^ ^ 1. This is because 
both the spin down and thermal time are less for the pure he- 
lium bursts (Azc cx 1 //i and ftherm oc l//i). As we discuss further 
in §5.2, Figure^ may explain the lack of burst oscillations seen 
in bursts with the characteristics of mixed H/He ignition. For 
these mixed H/He ignition bursts, ftherm Ai^ « 1 and any oscilla- 
tions in the flux from deeper regions can get washed out. 

3.5. Radius Expansion Bursts 

Many of the bursts which have shown oscillations are radius 
expansion bursts (Table in which super-Eddington luminosi- 
ties result in expansion of the photosphere to radii ~ 20-100 km 
(for a review see Lewin, van Paradijs, & Taam 1993). Burst 
oscillations have not been observed during the peak of radius 
expansion bursts, but are sometimes seen during the burst rise, 
and often during the tail of the burst, once the photosphere has 
fallen back to the stellar radius R (for an example, see Figure 
2 of Strohmayer et al. 1998c which shows a radius expansion 
burst from 4U 1636-54). 

We did not model radius expansion in §||. However, several 
of our models have super-Eddington fluxes, and thus are appro- 
priate for the early stages of these bursts. Because the opacity 
deep in the atmosphere is smaller than at the photosphere, the 
hydrostatic structure we calculate breaks down only in the very 
upper layers near the photosphere. Our radiative models are 
also appropriate for the cooling tail of these bursts, when the 
photosphere has fallen back to the stellar surface. 

Whatever the cause of the asymmetry on the neutron star sur- 
face, the fact that oscillations are not seen during the peak of 
radius expansion bursts is as we would expect, given the dis- 
cussion of §3.4. The long thermal time across the extended en- 
velope during radius expansion (Paczynski & Anderson 1986), 
as well as the fact that the horizontal and vertical lengthscales 
are similar, will hide any asymmetry and wash out the oscilla- 
tion. 

4. ANGULAR MOMENTUM TRANSPORT AND RECOUPLING 

We have shown that if different layers of the atmosphere con- 
serve their angular momentum, expansion during the burst re- 
sults in differential rotation within the burning layers. Thus as 



we noted in §3.4, the fact that a single coherent frequency is 
observed imp ies that either the burning layers are geometri- 
cally thin, or that they must rotate rigidly. While convection 
may enforce rigid rotation in the early stages of a burst, we do 
not e xpec t the atmosphere to be convective in the cooling tail 
(see §2.2). In this section, we investigate mechanisms which 



might transport angular momentum between the differentially- 
rotating burning layers. In addition, we investigate whether the 
burning layers can remain decoupled from the underlying cold 
ashes for the « 10 s duration of the burst and the implied several 
phase wraps. 

4. 1 . Kelvin-Helmholtz Instability 

Shear layers are notoriously unstable to hydrodynamic insta- 
bihties, in particular the Kelvin-Helmholtz instability (Chan- 
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drasekhar 1961). The shear may be stabilized by buoyancy, 
however, if the work that must be done against gravity to mix up 
the fluid is greater than the kinetic energy in the shear The im- 
portance of buoyancy is measured by the Richardson number, 
'Ri = / {dU / dzf (for example, see Fujimoto 1988), where 

is the Brunt-Vaisala frequency, a measure of the buoyancy. 
The Kelvin-Helmholtz instability occurs when Ri < 1 /4 (Chan- 
drasekhar 1961; Fujimoto 1988). 

There are two sources of buoyancy in the atmosphere, ther- 
mal buoyancy and buoyancy due to the composition gradient. 
We write the Brunt-Vaisala frequency as 



dlny ) 



Xp dz Xp dz 



(24) 



(Bildsten & Gumming 1998), where xq = d\nP/d\nQ with 
the other independent thermodynamic variables held constant. 
The first term of equation (Q) is the thermal buoyancy. For 
an ideal gas atmosphere carrying a constant heat flux in which 
the opacity is Thomson electron scattering opacity, this is A^^, = 
3nmpg^/2QkBT (Bildsten 1998), giving 



N,h ^3x 10"^ Hz f ^1 \ 1/2 

2^ ~ 7^1/2 



We estimate the composition gradient terms as A^^ 
gAln/x/i/, giving 



^«7xlO^Hz(AlnM)'/^(l^ 



1/2 



(26) 



where we have taken the density scale height H as the length- 
scale over which the mean molecular weight changes by 
Aln^-O(l). 

We show the Brunt-Vaisala frequency during the burst for 
two fully-radiative models in the top panel of Figure ^ We 
show the mixed H/He ignition model with F = 0.83 pEdd, and 
the pure He ignition model with F = 0.59 pEdd (see Table Q). 
We show the total buoyancy by a solid line, the thermal buoy- 
ancy as a dotted line and the composition piece as a dashed 
line. For the pure He ignition model, there is a peak in the 
buoyancy at the place where the hydrogen runs out (y w yd), at 
this depth the composition piece of the buoyancy dominates the 
thermal piece. Our models do not include the composition jump 
at the base between the burning layers and the ashes. The point 
marked by a cross at the base of each model shows an estimate 
of the buoyancy at the base, where we estimate the composition 
piece of the buoyancy using equation ( p6[ ) and take the ashes to 
have /i = 2.1 (for a single species, /i = A/l +Z). 

The Richardson number Ri = /{dU /dz)^ is plotted in the 
lower panel of Figure ^ for i^s = 300-600 Hz. Since Af] « 
i2z/R)n, we estimate dU/dz = Rd An /dz « 2n, giving 



Ri « 7000 



V300 Hz J 



N 



5 X 104 Hz 



(27) 



which agrees well with Figure ^. As suggested by Bildsten 
(1998), the strong buoyancy in the atmosphere gives Ri 3> 1. 
Thus we do not expect the Kelvin-Helmholtz instability to op- 
erate either within the burning layers or between the burning 
layer and the underlying ashes. 



4.2. Ekman Pumping 

We now investigate how fast viscosity acts to smooth out dif- 
ferential rotation. The molecular viscosity in the atmosphere is 
determined by ion-ion collisions, giving 



(28) 



(Spitzer 1962), where Ty = T/IO'^K, ps = p/lO^gcm-^ and 
log A is the Coulomb logarithm. Thus the time for molecular 
viscosity to transport angular momentum over a scale height by 
diffusion is 




«3.2h 



H 
10 m 



(29) 



much longer than a rotation period. In this case, it is possible 
to exchange angular momentum on a faster timescale than 
by the process of Ekman pumping. This mechanism involves 
a secondary circulation in which fluid elements are exchanged 
between the bulk of the fluid and the thin viscous boundary 
layer Ekman pumping is well-studied in fluid dynamics (Ben- 
ton & Clark 1974), and Livio & Truran (1987) suggested that 
it may operate in accreting white dwarfs. For a non-stratified 
fluid, the Ekman spin up or spin down time is ?£ « (t^/fiy^^ or 



f£ «2.5 s 



300 Hz 



1/2 



(30) 



(Greenspan & Howard 1963; Benton & Clark 1974). The 
thickness of the viscous boundary layer is « (v/fiy^^ « 
0.2 cm(iy/88cm2s-i)i/2(iy,/300 Hz)-'/^. 

Spin up in a stratified fluid is different since the buoyancy 
may inhibit vertical motion of fluid elements, limiting the ex- 
tent of the secondary circulation. Spin up in a cylinder of radius 
a with g parallel to 57 was first studied by Holton (1965; see 
Benton & Clark 1974 for a review). In this case, the secondary 
circulation is confined to a layer of vertical thickness « a{Vl/N) 
(Walin 1969). Sakurai, Clark & Clark (1971) found a similar 
result for a spherical stratified fluid, for which Ekman spin up 
occurs in a thin layer of radial extent w R(il/N). The neutron 
star atmosphere during a burst has H < R(il/N) « 60 m. Thus, 
if the spherical result is applicable to a thin layer on the surface 
of a sphere, we would not expect the buoyancy to inhibit Ekman 
pumping. In addition, we have ftherm ~ Ie, so non-adiabatic ef- 
fects may be important, reducing the effect of the thermal buoy- 
ancy (Sakurai et al. 1971). Detailed calculations are needed to 
find whether Ekman pumping operates during the burst. 

4.3. Baroclinic Instability 

The baroclinic instability is a hydrodynamic instability in 
which the fluid motions are close to horizontal, so it occurs 
even when the Richardson number is large. It has been well- 
studied in geophysics (Pedlosky 1987) and in astrophysics be- 
cause of its possible role in angular momentum transport in 
stellar interiors (Knobloch & Spruit 1982; Tassoul & Tassoul 
1982; Spruit & Knobloch 1984) and accreting compact objects 
(Fujimoto 1988, 1993). We start by showing that in the pres- 
ence of a vertical shear, the surfaces of constant pressure and 
density are inclined with respect to each other. This represents 
a store of gravitational energy which the baroclinic instabil- 
ity can tap. We then present the results of a stability analy- 
sis of a plane-parallel model, for which we find that the strong 
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buoyancy limits the baroclinic instability to short vertical wave- 
lengths 200 cm). 

4.3. 1 . The Nature of the Baroclinic Instability 

We first study the misalignment of the constant density and 
pressure surfaces that arises in the presence of vertical shear on 
a rotating star This result is well-known on the Earth, where 
differential heating between the equator and pole gives rise to 
the "thermal wind" (Pedlosky 1987). To investigate recoupling 
of the ashes and burning layers, we adopt a simple "two layer" 
model, in which we include only the buoyancy associated with 
the interface between the ashes and the burning layers. Thus, 
we take the upper (lower) layer to have constant density p+ (p-). 
Both layers are in hydrostatic balance in the vertical direction. 
As elsewhere in this paper, we assume the burning front has 
spread over the whole surface. 

We work in the rotating frame, in which the lower fluid is 
stationary. We assume the upper layer is rigidly rotating with 
angular velocity -A57 (so that Aft is positive for a spun down 
shell). Since it is moving, the upper fluid feels a Coriolis force 
in the transverse direction. In a time « 17"' « 10"^ s, a trans- 
verse pressure gradient will be established which balances the 
Coriolis force (geostrophic balance). However, in the lower 
fluid there is no transverse pressure gradient (it feels no Cori- 
olis force). Since the pressure must be continuous across the 
boundary, it cannot be horizontal, but must slope. 

Figure [I^(a) illustrates this. We show a small section of the 
two layers at latitude 6 (6 is the angle from the pole) in the (r, 6) 
plane. Because the star is rotating, the equilibrium isobars are 
not spherical, but nearly (the stellar radius is larger at the equa- 
tor than the pole by « R^ffi/g w 100 m for typical parameters). 
We adopt coordinates z and y perpendicular and parallel to sur- 
faces of constant pressure (dotted line). The vertical component 
of the rotation vector is il cos 9, and the upper fluid is moving 
out of the page with velocity U = RAil sin 6. 

Now consider the pressure changes as we move from point 
A to point B on the boundary. First we move through the 
lower fluid along the dashed line, horizontally a distance dy, 
then vertically upwards a distance dz- The pressure change 
along this path is dP- = -p-gdz- Second, we move through 
the upper fluid along the dashed line. This time, there is a 
pressure change while moving horizontally also, giving dP+ = 
-p+gdz- p+2nu cos 9 dy. Demanding dP- = dP+, gives the 
slope 

dz in^RAn P+ . . . 

— = — — cos6'sm6', (31) 

dy g Si Ap 

where Ap = p_ - p+ > 0. If we take the equilibrium equipo- 
tentials to be spherical (dy w Rd9) and integrate, we find the 
change in height of the boundary is Zb(9) = -ZbCos26+ const., 
where Zb = (n^R^ /2g)(An/n){p+/ Ap) or 



Zb = 35 cm f ^ — I 

V 300 Hz/ 



An/n 

4 X 10-3 



ApJ' 



(32) 



where we have taken gu = 2 and R= 10 km. 

The sloping interface represents a store of gravitational po- 
tential energy. Assuming constant density, the displaced mass 
per unit area is pzb{9), and the gravitational energy per unit area 
required to displace the interface is pgzb{9)^ /2. Integrating over 
a spherical surface, we find the gravitational energy stored in 
the interface is AEg = {\6/A5){ATTR^p+gzl/2){Ap/ p+), or 



AE, = 6 X lO^'* er: 



Ap 



Zb 



35 cm 



10^ g cm 



gl4 

2 



(33) 



The amount of displaced mass is k, AnR^plb ~ 4 x 10'^ g, 
so Eg is only a few keV per nucleon, much less than the en- 
ergy produced by nuclear burning. Thus establishing the slop- 
ing interface poses no energetic obstacle. For comparison, 
the energy in the shear is pHbU^/2 per unit area, giving 
Eu « {%nR^ /'i)pHb{RAVlf /2 or 



3 X lO^"^ erg 



10^ g cm 



Hb 
10 m 



V300 Hz/ 



An/n 

4 X 10-3 



(34) 



when integrated over the surface. 

The baroclinic instability acts to release the gravitational po- 
tential energy stored in the misaligned pressure and density sur- 
faces. Consider the fluid displacements shown in Figure [lO|(b). 
If the fluid element is displaced to point A, it is heavier than 
its new surroundings and is pushed back by the buoyancy of 
the interface. However, if it is moved to point B, it falls in the 
gravitational field, releasing energy. Any displacement within 
the so-called "wedge of instability" (Pedlosky 1987) is convec- 
tively unstable in this way. 

4.3.2. Results of Stability Analysis 

We have carried out a linear stability analysis of a plane par- 
allel two-layer model. We do not present our detailed calcula- 
tions, rather we summarize our results and use simple physical 
arguments to understand them. We start by considering motions 
about a latitude 9 = 6*0, and adopting a local cartesian coordi- 
nate system {x^y^z) where the transverse coordinate x (y) is in 
the (f) (9) direction. We include the effect of sphericity using the 
"beta-plane approximation" of geophysics (Pedlosky 1987), i.e. 
we write 

2r2 sin 61 w 217 sin 9o + {9- 9q) cos 6*0, (35) 

where we assume (9-9o) = y/R 1. We consider a chan- 
nel centered on 9() = tt/4 and of width ttR/4, stretching from 
the pole to the equator We write the fluid displacement as 
i£,x^^y^£,z)^ ^nd look for solutions cx exp{at + ikiX+ikyy) where 
kj, = 2Tr/Xx and ky = 27r/Ay are the transverse wavenumbers, 
and a is the growth rate. Pedlosky (1987) performs a similar 
analysis, although restricted to the case Ap/ p ^ 1. 

We find that small wavelengths are stable because the fluid 
displacements lie outside the wedge of instability. To see this, 
assume that the perturbations are in geostrophic balance, in 
which case the continuity equation gives £,z/X, « Ro^^/A_l, 
where Ro ^ U /2nXj_ « AO/20 is the Rossby number of 
the perturbation, and the transverse wavelength A^ = 2TT/k±, 
with k\ = k^. + ky. Then the angle of fluid displacement is 

&J£.y ~ Cc/C± ~ RAnH/X\^, where the vertical wavelength 
is set by the vertical extent of the burning layers A-. « H. 
For instability, the angle of fluid displacement should be less 
than the slope of the interface (eq. pi]]), requiring A^^ ^ 

(gHb/n^)iAp/p^), or 



A± > Abc = 7.5 X 10^ cm 



Hb 
10m 



1/2 



1/2 



/Ap 
V P+ 



1/2 



300 Hz. 



(36) 
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where we insert the correct prefactor Thus the nearly hori- 
zontal unstable displacements require transverse wavelengths 
of order the stellar radius or greater 

However, very large transverse wavelength perturbations can 
be stabilized because the vertical component of fl changes sig- 
nificantly across a wavelength. The changing vorticity provides 
a restoring force in the 6 direction (Pedlosky 1987). This is 
the same force which supports Rossby waves (Pedlosky 1987; 
Brekhovskikh & Goncharov 1994; Dutton 1995). We find that 
the convective instability overcomes the Rossby wave restoring 
force when the shear velocity is greater than the Rossby wave 
speed, U > Uro ~ 2nX\/R (see also Pedlosky 1987). Thus we 
require A_l ^ R(An/ny/^, or 

A,<A.„ = 3.3xlO-^cm(^)'''(^) (37) 

for instability, where we insert the correct prefactor 

Transverse wavelengths which satisfy equation ( |3^ do not 
satisfy equation (^7|). Thus in the context of our plane-parallel 
model, we conclude that, when H = 10 m, large scale mixing 
between the burning layers and ashes by the baroclinic insta- 
bility does not occur Instability can occur when Abc < ^Ro, or 
Ap/p+ < {4n^R/g)(An/n){R/H). since AQ/n w 2H/R, this 

so mixing may occur as the burning layers cool and the den- 
sity contrast with the ashes decreases. The growth rate of 
an unstable mode is oc k^-, so that writing ky- « m/R (so that 
e'**' ^ e"""^), we estimate the fastest growing modes will have 
m « 20 (A^/3 X 10^ cm). Several authors have studied two 
layer models on a sphere (Hollings worth 1975; Hollings worth, 
Simmons & Hoskins 1976; Simmons & Hoskins 1976, 1977; 
Moura & Stone 1976; Warn 1976), but for parameters of inter- 
est for the Earth, namely Ap/p < 1 and Ail/Q « 0.01. They 
find good agreement with growth rates calculated with beta- 
plane models. However, further calculations are needed to ex- 
tend our analysis to the spherical case, and determine the spher- 
ical eigenfunctions and growth rates. 

4.3.3. Short Wavelength Modes 

So far we have considered only vertical wavelengths A, « H, 
because of the simple vertical structure of the two layer model. 
However, if we allow perturbations with A. < H, we expect 
to find instability for small enough vertical wavelengths. This 
is because for small vertical wavelengths, the fluid displace- 
ment is able to lie within the wedge of instability (Figure |l^ 
while still having a transverse wavelength which is unaffected 
by Rossy wave restoring forces (eq. [p7[]). 

We make a simple estimate of the vertical wavelength at 
which instability occurs by repeating the arguments leading to 
equation (p6|), but this time allowing A, < Hh- In this case, Xbc 
is given by equation ( |3^ ) with H/, replaced by A,. For instability, 
we require Abc < Aro (eq. [|7[]) giving 

A-<200cmf^V^)f^V (39) 
V^x 10-V ^300Hz/ V^p/ 

where we take = 2 and R= IQ km. We expect modes with 
vertical wavelengths small enough to satisfy equation (^9|) to be 
unstable. 



Fujimoto (1988, 1993) used the short wavelength baroclinic 
modes to define a turbulent viscosity i^turb A?/ fgrow where the 
growth time of the instability is fgjow ~ v'^Ri^^^, or 




(Pedlosky 1987; Fujimoto 1988). The time to transport an- 
gular momentum across a scale height is « H^/i'imh ~ 
(H/X^)^, or 

t^^6s(^)(JL.V(^l^Y (41) 
Vo.25s/VlOmy V200cmy ^ ' 

Thus turbulent transport of angular momentum driven by the 
short wavelength baroclinic modes could be important during a 
burst. In particular, these could act to force the burning shell to 
rigidly rotate. 

4.4. Magnetic Field Winding 

The majority of neutron stars in LMXBs show no evidence 
for coherent pulsations in their persistent emission. This im- 
plies that these neutron stars are weakly magnetic (B <C 10^ G) 
so the accretion flow is likely not disrupted before it reaches 
the neutron star surface. We now ask what effect would a weak 
magnetic field have on the shearing atmosphere during a burst? 
The Ohmic diffusion time across a scale height is t = 4TTaH^ /c^, 
or 

lO^s f H Yf T 8 \ 

fd.ffuse«— (^Y^j (^y^j ^42^ 

(Brown & Bildsten 1998) where we estimate the conductivity 

as 

^ 2(2kBT)y^ , ... 

TT^I'^mJ Ze^logA 

(Spitzer 1962), Z is the mean ionic charge and log A is the 
Coulomb logarithm. Since fdiffuse 3> 10 s, the MHD limit ap- 
plies during the burst and the shearing atmosphere will bend 
the magnetic field lines. An initially poloidal magnetic field 
would prevent the shearing if the energy density in the field is 
greater than the shear energy, B^/87r > p(RAil)^ /2, or 

B>10-Gpr(^)(^). (44) 

This is a similar limit to those from lack of persistent pulsations 
and from spectral modelling (Psaltis & Lamb 1998). The ob- 
served shearing for ten seconds most likely rules out fields this 
strong. 

What can we say about a weaker, initially poloidal field like 
that often discussed for the the progenitors of millisecond pul- 
sars, Z? 10** - 10^ G? Fields this weak can also have a sig- 
nificant effect when in the ideal MHD limit. In this case, the 
horizontal shearing and dragging of the field lines generates a 
strong toroidal field. Spruit (1999) recently considered the in- 
teraction of magnetic fields and differential rotation in stellar 
interiors. He considered an initially poloidal field (B.,Bg,Q), 
and showed that grows with time according to = NB^ 
where is the number of windings, N = RsinO J \'Vil\dt. We 
have Vn « An/H, giving « {RAQ/H)t « 2^t. Thus we 
find B^ K, B, after only a few revolutions of the star, much less 
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than the wrap around time of the shear (« 1 s). The toroidal 
field grows so quickly because the field is sheared on a ver- 
tical scale H ^ R. Another way to see this is to consider a 
vertical tube of fluid in the shearing region (with height H and 
cross-section area A,) that encloses a magnetic flux $ = B^A,. 
After one differential wrapping time (roughly one second), 
this fluid element will be stretched into a very thin tube of 
length IttR and so will have shrunk in cross-sectional area to 
Af = AjH /2ttR. If flux-freezing holds, then after one second, 
B0 « BAi/A f w B,2ttR/H « lO^B,. 

This simple argument says that an initial field of w 10'' G can 
become dynamically important during the burst. However, fur- 
ther theoretical investigations are needed to see if the implied 
wrapping of the field is possible without instabilities setting in. 
If correct, this result implies that the bursters with frequency 
drift have magnetic fields much weaker than the B ^ 10**-10'^ G 
fields usually assumed. 

4.5 . Summary of Coupling Mechanisms 

We have investigated a number of different mechanisms that 
could enforce rigid rotation within the burning layers, or act to 
couple the burning layers to the underlying ashes. We find that 
the strong buoyancy of the atmosphere prohibits the Kelvin- 
Helmholtz instability and mixing between the ashes and the 
burning layers by the baroclinic instability. Short vertical wave- 
length baroclinic modes within the burning region may be un- 
stable and vigorous enough to force the burning layers to rigidly 
rotate. The timescale for viscosity to act via an Ekman pump- 
ing mechanism may be a few seconds or less (eq. [Ml), but 
the effect of buoyancy, in particular from the composition gra- 
dients, is not clear. Finally, our initial estimates suggest that 
winding of a weak magnetic field by the differential rotation is 
important during the burst. 

In summary, we have not found a robust hydrodynamical 
mechanism that recouples the burning layers to the star during 
the R:! 10 s burst duration. It thus seems plausible that the burn- 
ing layers remain decoupled throughout the burst. Short wave- 
length baroclinic instabilities may operate within the burning 
layers, although more work is needed to see whether they cause 
the burning shell to rigidly rotate. 

5. DISCUSSION 

We now summarize the hydrostatic expansion of the atmo- 
sphere during a Type I X-ray burst, and the implied spin evo- 
lution of the decoupled burning layers. We find the expected 
spin changes are of the same order as the observed frequency 
drifts, thus supporting the picture proposed by Strohmayer et al. 
(1997a). We then use our results to address two observational 
issues: (i) why oscillations are not seen in all bursts, and (ii) the 
observed long term stability of the oscillation frequency. 

5.1. Temporal Evolution of the Burning Shell 

The atmospheric evolution during the burst is shown 
schematically in Figure [ll] for (a) mixed H/He ignition and (b) 
pure He ignition. As revealed by one-dimensional simulations, 
the rapid energy release from helium burning reactions makes 
the initial stages of most bursts convective. Later in the burst, 
the convection zone retreats, leaving a purely radiative atmo- 
sphere. Figure [l^ shows the thickness of the atmosphere which 
contains 90% of the mass, Az(90%), as a function of the flux for 
constant flux radiative atmospheres (solid and dotted lines) and 
for the convective models of Table || (solid and open squares). 



For the radiative models, the upper curve is for a composition 
profile the same as at ignition, while the lower curve is for a 
composition of ^■'Kr (^^Ni) for the mixed H/He (pure He) case, 
chosen to represent the products of burning. Initially, the evo- 
lution of the thickness of the layer is along or above the up- 
per curve as it ignites and heats up, depending on the extent of 
the convection zone. As nuclear burning proceeds, the mean 
molecular weight increases, and the thickness decreases, even- 
tually moving back along the lower curve as the atmosphere 
cools. Just before ignition (upper curve, F <C /^Edd), the thick- 
ness of the accumulated layer is « 5 m. During the burst, the 
atmosphere expands hydrostatically by Az w 10^0 m (« 5- 
20 m) for mixed H/He (pure He) ignition. For a given flux, the 
convective models have a greater thickness than fully-radiative 
models because the temperature profile is steeper in the con- 
vection zone, giving a larger base temperature (compare Tables 
[3] and H). As we showed in §2.2, radiation pressure plays an 
important role when F « F^di, acting to lower the density and 
increase the hydrostatic expansion. 

We showed in §^ that in a radiative atmosphere, the mag- 
nitude of the spin down depends on depth. Thus to ob- 
serve a single, coherent frequency requires either the atmo- 
sphere to be rigidly rotating, or the burning region to be ver- 
tically thin (Miller 2000), so the differential rotation across it 
is small. Convection may enforce rigid rotation; however, one- 
dimensional calculations of bursts show that convection persists 
only for the first <; 1 s of the burst when the energy production 
is dominated by temperature sensitive and rapid helium burn- 
ing reactions. In §4, we found that short wavelength baroclinic 
instabilities may act to enforce rigid rotation within the burning 
layers. Theoretical studies of detailed burst models and further 
investigations of angular momentum transport mechanisms are 
needed to determine whether a coherent pulse is possible once 
convection has ceased. 

To compare the predicted spin down with observed values, 
we assume that some mechanism operates to enforce rigid ro- 
tation within the burning layers, but that the burning layers re- 
main decoupled from the underlying ashes, as suggested by the 
results of §4. Figure |l^ shows the spin evolution of the at- 
mosphere assuming that the whole atmosphere rigidly rotates. 
For both convective and fully-radiative models, we ignore any 
differential rotation, and assume IVt is constant, where / is the 
moment of inertia of the atmosphere (eq. [p^). As in Figure 



12, the evolution during the burst is shown by the arrows, first 



the atmosphere spins down as flux increases, then spins back 
up as burning proceeds, increasing the mean molecular weight, 
and as the atmosphere cools. To compare with observations, 
we also indicate the observed frequency shifts (Table |l]). For 
those bursts in which the oscillation frequency was seen only in 
the tail (Table |l]), we plot AO/0 as a lower limit, since some 
cooling will have occurred before the oscillation is first seen. 

Figure |l3] shows that the change in spin frequency for pure 
He ignition models is less than that for mixed H/He ignition, 
because of the smaller hydrostatic expansion. For the largest 
frequency shifts observed, AO/O « 8 x 10"^, pure He ignition 
models require a base temperature very close to the limiting 
value from radiation pressure (eq. [|l^]) and F > F^dd, so that 
the density at the base is decreased, enhancing the expansion. It 
is much easier to get Aft / ft values this large with mixed H/He 
ignition, although as we describe below, all bursts with oscilla- 
tions so far have the characteristics of pure He bursts. 

In §111 we found that the time for heat to diffuse from the 
burning layers to the photosphere is « 1 s. This delay explains 
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why the spin down of the atmosphere is not observed at the 
beginning of a burst; the expansion occurs before the signal is 
seen. We have also shown that the vertical propagation of heat 
through the atmosphere can affect the oscillation amplitude. A 
requirement for a large amplitude oscillation is that the time to 
transport heat from the burning layers (or the top of the con- 
vection zone) to the photosphere must be small compared to 
the shearing time across the atmosphere « Aj/"'. If so, the 
shearing atmosphere transmits the burning pattern from below, 
much like passing a flashlight under a piece of paper. If not, the 
oscillation will be smeared out, the atmosphere being heated 
uniformly by the quickly revolving burning layers. 

Our results show that there are important differences in spin 
evolution depending on whether the helium ignites in a pure he- 
lium or mixed H/He environment. Because of the greater mean 
molecular weight of pure helium as opposed to a solar com- 
position, the expansion of the atmosphere in pure He bursts is 
roughly half that of mixed H/He bursts, giving a smaller Afl/il 
(Figure [l^). Also, the ratio of thermal time to shearing time 
is smaller for a pure helium atmosphere, and the energy pro- 
duction more localized because of the temperature sensitivity 
of He burning reactions, increasing the likelihood of observing 
coherent oscillations for pure He bursts. 

5.2. Why Do Some Bursts Show Oscillations, But Not All? 

A puzzle is that oscillations are not seen in bursts from all 
sources, or in all bursts from a particular source. So far, there 
has been only one study of the relation between burst oscilla- 
tions and properties of the bursts. Muno et al. (2000) studied 
nine X-ray bursts from KS 1731-26 which occurred during ob- 
servations with RXTE at various times between 1996 July and 
1999 February. They found that bursts which show coherent 
oscillations are of short duration (« 10 s), show radius expan- 
sion, and have high peak flux. These characteristics are typical 
of helium rich bursts. Bursts with longer durations typical of H 
burning showed no oscillations. More studies are needed to see 
if this is a general result for other sources. Certainly all reported 
burst oscillations that we can find appear to be in bursts of short 
duration typical of pure He ignition. As far as we are aware, 
oscillations have not been detected in bursts of long duration 
(> 50 s) typical of H burning. 

Why should only pure-He ignition bursts show oscillations? 
The differences we find between pure He and mixed H/He 
bursts may provide an explanation. As we described in §3.4 
(see Figure @), we find that the thermal time is less and the 
shearing time is greater for pure helium bursts, making it more 
likely that a large amplitude oscillation from deep regions may 
propagate outwards. 

5.3. Long Term Stability of the Oscillation Frequency 

The frequency evolution reported for six bursts is well fit by 
the model v{t) = vq - Ave~'^'^ , where I'o is the frequency in the 
burst tail, Aiy the amount by which the oscillation frequency 
changes and r the decay time. The parameters of these fits are 
given in Table 1 . The observed decay times r « 2-7 s are simi- 
lar to the expected cooling time of the burning layers (eq. |]23|]). 
However, we found in §4 that there are mechanisms that might 
couple the burning layers to the star on a similar timescale. Can 
we distin guis h between these two possibilities observationally? 

Figure |l2| shows that if the burning layers remain decou- 
pled, the frequency observed in the burst tail will not be that 
of the neutron star spin. Because of the greater mean molecular 
weight of the ashes, the thickness of the cooling atmosphere in 



the tail is different to the initial thickness by « 1 m. This change 
in thickness depends on how complete the burning was, so that 
variations in the energetics and burning in the burst would trans- 
late into a scatter of « 1 part in 10'* in the final frequency. If 
the oscillation frequency could be shown to be more stable than 
this from one burst to the next, it would imply that the ashes 
and burning layers must recouple during the burst decay. This 
is complicated due to Doppler shifts from the orbital motion of 
the neutron star, which change the observed frequency at a level 
w 2 X 10~-'(voi-b/600 km s"') across the orbit. Strohmayer et al. 
(1998b) showed that two bursts from 4U 1636-54 showed a fre- 
quency difference consistent with orbital Doppler shifts, while 
two bursts from 4U 1728-34 separated by more than one year 
showed the same asymptotic frequency to « 1 part in 10"^ (see 
Table 1). If the orbital parameters are known and the orbital 
Doppler shift accounted for, more measurements such as these 
would indicate whether recoupling occurs in the burst tail. 

6. CONCLUSIONS 

We have shown that the hydrostatic expansion, spin down, 
and later spin up during contraction of the neutron star atmo- 
sphere naturally explains the magnitude and sign of the ob- 
served changes in the nearly coherent oscillations observed in 
Type I X-ray bursts. Our results support the simple picture pro- 
posed by Strohmayer et al. (1997a) to explain the frequency 
evolution seen in Type I bursts, and thus the identification of 
the burst oscillation frequency with the neutron star spin. 

The amplitude of the oscillation is set not only by the lat- 
eral extent of the asymmetry on the neutron star surface, but 
also by the vertical propagation of heat through the shearing 
atmosphere. In addition, we find that the spin behavior differs 
depending on whether the burst results from pure He ignition or 
He ignition in a H-rich environment. The spin down is smaller 
by roughly a factor of two for pure He ignition as opposed to 
mixed H/He ignition, because the pure He models expand out- 
wards less due to the larger mean molecular weight. In addition, 
the time for heat propagation through the atmosphere is smaller 
for pure He models. 

This might explain why oscillations have not been detected 
during bursts of long duration (^ 30 s) typical of H burning (for 
example, the bursts from GS 1826-24 shown in Figures 2 and 
5 of Kong et al. 2000). In particular, the recent study of Muno 
et al. (2000) of bursts from KS 1731-26 found oscillations only 
during those bursts with short durations « 10 s, typical of pure 
He ignition. However, the largest values of Ail/il that are ob- 
served (for 4U 1702-43 and 4U 1728-34) require base tempera- 
tures very close to the limiting value from radiation pressure for 
pure He bursts, and this may prove to be a problem. More ob- 
servational studies of the phenomenology of burst oscillations 
and their relation to burst properties and source properties, in 
particular accretion rate, are needed. With further theoretical 
work, these studies could give us important clues as to how the 
nature of nuclear burning during bursts depends on accretion 
rate, and help to resolve the disagreement with one-dimensional 
models first found with EXOS AT (Bildsten 2000 and references 
therein). 

We have not found a robust hydrodynamical mechanism 
that recouples the burning layers to the underlying star in the 
« 10 s burst duration. However, the coherence of the observed 
burst oscillations suggests that little differential rotation occurs 
within the burning layer itself. As we note in §Q this might 
be accomplished through short vertical wavelength baroclinic 
instabilities, though more theoretical work is needed. 
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We have not discussed the complex question of the lateral 
spreading of the burning front during the burst rise, or the 
cause of asymmetry in the cooling tail of bursts. In the sim- 
plest picture of a spreading burning front (Strohmayer et al. 
1997b; Strohmayer et al. 1998c) it is a mystery why oscilla- 
tions are seen in the burst tail, when the whole surface has ig- 
nited. One possible explanation is that an asymmetry arises 
because the atmosphere takes a finite time to cool once ig- 
nited, and a finite time to spread around the star. Thus the 
matter on the opposite side of the star from the ignition point 
will be a little hotter since it ignited later If the flux at each 
point is cx e"'/^ once ignited, and the burning front propagates 
around the whole star (starting at the equator) in a time fbum, 
it is straightforward to show that the peak-to-peak luminosity 
amplitude in the tail for an observer looking down on the igni- 
tion point is AL/L = 2(1 -exp(-fbum/T))/(l +exp(-fbun,/2T))^, 
or, for fburn ^ T, AL/L fburn/2T. Thus for a spreading time 
fburn ~ 0.5 s and r « 10 s, the peak-to-peak amplitude in the 
tail is a few percent, a little less than the 5-20% amplitudes 
observed (Smith, Morgan & Bradt 1997; Zhang et al. 1998; 
Strohmayer et al. 1998a). Another possibility is that the oscil- 
lations come and go because of changes in the vertical transport 
of heat rather than changes in the size of the asymmetry on the 
stellar surface. For example. Figure |^ shows that as the atmo- 
sphere heats up during the burst rise, the ratio of thermal time to 
wrap-around time ftheimAt^ increases. As this ratio increases, it 
becomes more likely that the oscillation will be washed out. It 
may be that at the burst peak, the oscillations disappear because 
of this effect. 

We have not addressed what mechanism might cause fre- 
quency doubling during bursts, as seems to be the case for at 
least one source (4U 1636-56; Miller 1999). There are two pro- 
cesses that we have discussed in this paper that may play a role. 
The first is the recoupling mechanism. If the burning layers 
recouple to the neutron star during the burst, they may do so 
in an m = 2 pattern, leading to frequency doubling. Our lin- 
ear instability analysis of the two-layer beta-plane model (§ 4-. 3 ) 
suggests that the baroclinic instability would have m w 20 when 
it goes unstable. However, we have not calculated the correct 
spherical eigenfunctions. The second possible process is the 
azimuthal wrapping of the the magnetic field. If an initially 
dipolar field is wound up, an m = 2 pattern might emerge in the 
flow. 

A test for the rotational modulation hypothesis might be the 
distribution of observed Afi/fi values. For sources undergo- 
ing X-ray bursts of the same type (for example, similar ac- 
cretion rates and same type of ignition) on similar mass neu- 
tron stars, we would expect the distribution of Ail/H. values 
to be the same for each source, and independent of the spin 
frequency or number of hotspots on the surface. In fact, the 
distribution of Ail/il so far (Table 1) looks bimodal, though 
we should be wary of small number statistics. Two sources, 4U 
1728 and 4U 1702, which have i/q « 350 Hz show frequency 
changes Aiy/iy w 5-8 x lO""*, whereas the other four objects 



which have j/q ~ 550 Hz show frequency changes about half 
as large, Av jv w 1-4 x lO"-'. If further observations show this 
bimodal distribution to be true, this might indicate that the ex- 
pansion of the atmosphere is less for the i/q w 550 Hz objects, 
for example, because of pure He rather then mixed H/He igni- 
tion. Alternatively, we note that the Av values are very similar 
for all six objects, and this might point to another explanation 
for the observed frequency evolution. Care must be taken that 
the observed Af2 /il represents the change in burst oscillation 
frequency from the burst rise to the burst tail. For example, if 
the oscillation is seen only during the burst tail, the observed 
A57/f2 does not represent the full spin down of the atmosphere, 
as in this case some cooling will have occurred before the os- 
cillation is seen. 

We have not discussed other explanations for the frequency 
drifts. One possibility is that the burst oscillation is a non- 
radial oscillation (NRO) on the neutron star surface (Livio & 
Bath 1982; McDermott & Taam 1987; Bildsten & Cutler 1995; 
Strohmayer & Lee 1996). Non-radial oscillations of the fre- 
quency observed would be plausible on a rapidly rotating star, 
as the Coriolis force takes the typically much lower-frequency 
modes and makes them mugh higher (Bildsten, Ushomirsky 
and Cutler 1996). The observed frequency would be close to 
the spin frequency of the star or some multiple mfi, and might 
change by a 1-2 Hz during the burst as the atmosphere cooled. 
However, it is not clear why only a single mode would be se- 
lected, and why that frequency would be stable on timescales 
of a year or more, for which conditions in the surface layers of 
the star would most likely change significantly. 

The physics of expansion-induced shear during shell flashes 
may prove important in other contexts. For example, mixing 
due to shear instabilities during a classical nova explosion could 
provide the enrichment of CNO elements necessary to achieve 
a fast nova (see Truran 1982 for a review of the enrichment 
problem). Previous authors have investigated shear mixing dur- 
ing accumulation of the fuel (see Livio & Truran 1990 for a 
review), but not during the burst itself when shear is induced by 
expansion of the atmosphere as in a Type 1 X-ray burst. A sim- 
ilar situation arises in AGB stars during He shell flashes, where 
mixing of carbon into the proton rich envelope is necessary to 
allow production of s-process elements (Iben 1991; Sackmann 
& Boothroyd 1991; Langer et al. 1999), and perhaps could be 
achieved by shear instabilities during shell flashes. 
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Fig. 1 . — Profiles of the atmosphere immediately prior to (a) mixed hydrogen and helium ignition and (b) pure helium ignition. The hatched region shows where 
helium ignition occurs for helium mass fraction Y between 1.0 and 0.3. All models have a base flux of Fi, = 150 keV per nucleon (see text), (a) The solid lines show 
the profiles for Zcno = 01 and (bottom to top) ni/mEdd = 0.03,0.1,0.3. The dashed lines show the m = 0.1 niEdd case for Zcno = 0.005 (bottom curve) and 0.02 
(top curve), (b) The solid hnes show the profiles for Zcno = 0.01 and (bottom to top) Jfi/"'Edd = 0.01,0.015 J). 02. The dashed lines show the m = 0.01 niEdd case for 
Zcno = 0.005 (top curve) and 0.02 (bottom curve). The black dots show where the hydrogen runs out (eq. M). 
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Fig. 2. — Temperature profiles (upper panel) and spin down (lower panel) for the convective models of Table g. In each case, the black dots mark the top of the 
convection zone at column depth yc (Table pi), (a) We show models for mixed H/He ignition with m = 0.1 njEdd, Zcno = 0.01, and V6 = 2.04 X 10* gcm"^. These 
models have {yi,-yc)lyb = 0.5,0.8,0.95 andTD.99, Ti/lO' K= 1.2, 1.5, 1.7 and 1.7, and F/fEdd = 0.31.0.74, 1.15 and 0.89. (b) We show models for pure He ignition 
with m = 0.015mEdd, ZcNO = 0.01 and y,, = 2.53 X lO^gcm"^ -p^ese models have (yb-yc)/yb = 0.5,0.8,0.95 and 0.99, rfc/lO'K = 1.3, 1.6, 1.95 and 2.0, and 
F/Fem = 0.35,0.59, 1.21 and 1.06. The difference in slope of the convective adiabat between models is due to the differing contribution of radiation pressure. 
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Fig. 3. — Density as a function of height (solid lines) for the convective models of Table 3 and Figure ^for (a) mixed H/He ignition and (b) pure He ignition. The 
density profile just before ignition is shown as a dashed line in each case. The black dots mark the height which encloses 90% of the mass. In (a) the density profiles 
of the two models with f /Fedd = 0-^9 and 1.15 overlap. 
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Fig. 4. — Temperature profiles (upper panel) and spin down (lower panel) for the fully-radiative models of TableW. (a) We show constant flux models for mixed 
H/He ignition with m = 0.1mEdd, Zcno = 0.01, andyj, =2.04 X lO^gcm"^ ^hese models have F/^Edd = O.SLO.V'Tand 1.15 and Ti/lO' K= 1.12,1.35 and 1.49. 
(b) We show constant flux models for pure He ignition with m = 0.015 mnddi Zcno = 0.01, and yi, = 2.53 X 10* gcm"^. These models have F/FEdd = 0.35,0.59 and 
1.21 and Ti/lO'' K= 1.19, 1.34 and 1.57. 




Fig. 5. — Density as a function of height (solid lines) for the fully-radiative models of Table 4 and Figure ^for (a) mixed H/He ignition and (b) pure He ignition. 
The density profile just before ignition is shown as a dashed line in each case. The black dots mark the height which encloses 90% of the mass. 
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Fig. 6. — The AQ/Q profiles at different times as the convection zone shrinks, for mixed H/He ignition with m = O.lmEdd and Zcno = 0.01. We take the initial 
extent of the convection zone to be (yb-yc)/yb = 0.8 and the base temperature ri,/10'K= 1.5. As the convection zone retreats, we assume the base temperature 
decreases linearly (eq. |p2|) to T^/ lO' K = 1 .0. We show Af2 as a function of depth for the initial model (80% of mass convective), an intermediate case (50%) 
and when the convective zone has almost vanished (5%). These models have 7;,/10' K = 1.5, 1.31, and 1.03 and F/F^^i = 0.74,0.49, and 0.21 respectively. 
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Fig . 7 . — Models of the cooling atmosphere in the tail of the burst. We show the physical thickness which encompasses 90% of the mass Az(90%), as a function 
of the flux F/Fem for constant flux atmospheres. The solid curve is for the ni = 0.1 fJiEdd, Zcno = 0.01 mixed H/He ignition model, where we take the composition 
to be pure '^Kr. The dashed lines are for pure helium ignitions, m = 0.01 mndd, Zcno = 0.01 (upper curve) and m = 0.015 niEdd, Zcno = 0.01 (lower curve). In this 
case, we take the composition to be '^Ni. In each case, we show Az(90%) just before ignition by a horizontal solid line. 
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Fig. 8. — Thermal time ?therm (eq.[p3||) and wrap around time i^v ' for tlie convective models of Table H. In each case, the solid line shows ^themn and the hatched 
region Ai/"' for fs=30O-600 Hz. (ajMixed H/He ignition convective models with (top to bottom) f /fEdd = 1.15,0.74 and 0.31. (b) Pure He ignition convective 
models with (top to bottom) f /^Edd = 1.21,0.59 and 0.35. 
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Fig. 9. — Brunt- Vaisala frequency N /Itt and Richardson number Ri for (a) the fully-radiative, constant F = 0.74 Fem mixed H/He ignition model and (b) the 
fully-radiative, constant F = 0.59 F^jd pure He ignition model. We show the total buoyancy by a solid line, the thermal buoyancy as a dotted line and the composition 
piece as a dashed line. The cross shows an estimate of N/2tt at the base, N^^^ = (g / Hi,)(Aln n), where we take the ashes as having fi = A/(l +Z) = 2.1. For the 
pure He ignition model (b), there is a peak in the buoyancy at the place where the hydrogen runs out (y « yd), at this depth the composition piece of the buoyancy 
dominates the thermal piece. 




< pole equator > 

Fig. 10. — A section of the two-layer model. The burning layers and ashes are represented by layers of different constant density. The radial coordinate z increases 
upwards, the transverse co-ordinate y increases to the right. The vertical component of the rotation vector is Q cos 6, and the upper fluid (density p+ < p-) is moving 
out of the page with velocity RAQsin(9. (a) The boundary between the two fluids (solid line) slopes because in the upper fluid there is horizontal pressure gradient 
which balances the Coriolis force (see text), (b) The nature of the baroclinic instability. If the fluid element is moved to point A, it is heavier than its surroundings, 
and experiences a restoring force. If it is moved to point B, however, it is heavier than its new surroundings but has dropped in the gravitational field, releasing 
energy. AH motions within the so-called "wedge of instability" are convectively unstable in this way. 
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Fig. 11. — Schematic of the evolution of the atmosphere during a burst for (a) mixed H/He ignition and (b) pure He ignition. The rapid energy release from 
helium burning reactions makes the initial stages of many bursts convective. After a time < 1 s (see text for references to simulations of bursts), the convection zone 
shrinks and disappears, leaving a radiative atmosphere. For pure He ignition, the burst duration is typically 10 s; for mixed H/He ignition, rp process H burning 
in the tail can prolong the burst for a further > 30 s. 
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Fig. 12. — The thickness of the atmosphere which ccmtains 90% of the mass, Az(90%), as a function of the flux for constant flux radiative atmospheres (solid 
and dotted lines) and for the convective models of Table bI (solid and open squares). For the radiative models, the upper curve is for composition profile the same as 
at ignition, while the lower curve is for a composition of^^Kr ('*Ni) for the mixed H/He (pure He) case, chosen to represent the products of burning. Initially, the 
evolution of the thickness of the layer is along or above the upper curve as it ignites and heats up, depending on the extent of the convection zone. As nuclear burning 
proceeds, the mean molecular weight increases, and the thickness decreases, eventually moving back along the lower curve as the atmosphere cools. Just before 
ignition (upper curve, F <g; FEdd)i the thickness of the accumulated layer is 5 m. During the burst, the atmosphere expands hydrostatically by Az 10^0 m 
(Ri 5-20 m). In the cooHng tail, the thickness may differ by 1 m from that before the burst. 
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Fig. 13. — The spin evolution of the atmosphere for convective (squares) and radiative (lines) models, assuming rigid rotation is maintained-lhroughout the 
atmosphere. We plot the fractional spin frequency change of the burning shell, Af2/f2, such that Af2/f2 > indicates spin down. As in Figure Il2l the evolution 
during the burst is as shown by the an'ows. We also indicate the observed frequency shifts for those objects listed in Table Im by horizontal bars. For those bursts in 
which the oscillation frequency was seen only in the tail, we plot the AQ/Q value as a lower hmit. 



